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ABSTRACT 

In this Letter, we investigate how galaxy mass assembly mode depends on stellar mass M„, using a large 
sample of ^10, 000 low redshift galaxies. Our galaxy sample is selected to have SDSS Rgo > 5".0, which 
allows the measures of both the integrated and the central NUV-r color indices. We find that: in the M,* - ( 
NUV-r) green valley, the M* <10*° Mq galaxies mostly have positive or flat color gradients, while most of 
the M* > 10'°'^ Mq galaxies have negative color gradients. When their central D„4000 index values exceed 
1.6, the M* < 10'°° Mq galaxies have moved to the UV red sequence, whereas a large fraction of the M* > 
I 010.5 galaxies still lie on the UV blue cloud or the green valley region. We conclude that the main galaxy 
assembly mode is transiting from "the outside-in" mode to "the inside-out" mode at M* < 10 '° Mq and at 
M* > 10 '° ^ Mq. We argue that the physical origin of this is the compromise between the internal and the 
external process that driving the star formation quenching in galaxies. These results can be checked with the 
upcoming large data produced by the on-going lES survey projects, such as CALIEA, MaNGA and SAMI in 
the near future. 

Subject headings: galaxies: evolution - galaxies: star formation 


1. INTRODUCTION 

The stellar mass of a galaxy, M*, is one of its most fun¬ 
damental properties. The assembly of M* is a long-standing 
issue in the field of galaxy formation and evolution. To elu¬ 
cidate how M* evolves, one would expect to recover the star 
formation history (SEH), both in space and time, for individ¬ 
ual galaxies. In practice, the SEH of a galaxy is determined 
by finding the most plausible combination of evolved single 
stellar populations (SSP) that matches its observed spectrum 
or spectral energy distribution (SED). This so-called "fossil 
record method" has been widely applied to the integral field 
spectroscopic (lES) data and the multi-band images with good 
spatial coverage to study galaxy assembly (Kong et al. 2000; 
Lin et al. 2013; Perez et al. 2013; Cid Eernandes et al. 2013). 

Previous works have found two distinct stellar mass assem¬ 
bly modes. They are termed as "the inside-out" and "the 
outside-in" modes in the literature (Sanchez-Blazquez et al. 
2007; Perez et al. 2013). In the "inside-out" scenario, mass 
assembly is firstly finished in the galactic central region. Eor 
instance, in a young disk galaxy, the disk instability will in¬ 
duce gas inflow and trigger starburst in the galactic center as it 
evolves. Besides of star formation, the inflow gas can also fuel 
the central supermassive blackhole. The subsequent AGN 
feedback (such as blowing the gas out of the galaxy, or heat¬ 
ing it against cooling down to form new stars) will then sup¬ 
press the central star formation, leaving a compact, quiescent 
galactic bulge (Dekel& Burkert 2014; Bournaud et al. 2014). 
Compared to the central part, its outskirts form through much 
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gentler process, such as the star formation driven by the grad¬ 
ually accretion of cold inter-galactic medium (IGM), or the 
accretion of small satellite galaxies. Many works found that 
massive spiral galaxies have negative color or stellar age gra¬ 
dients, supporting this picture (Wang et al. 2011; Perez et al. 
2013; Lin et al. 2013; Sanchez-Blazquez et al. 2014). Com¬ 
pared to their massive counterparts, low mass galaxies have 
much shallower potential wells, making them less capable to 
retrieve the gas blowing by the galactic outflows, or accrete 
cold gas to form new stars. In addition, the role of AGN 
feedback may not be so efficient in low mass galaxies, as 
their AGN occupation fraction is very low (Kauffmann et al. 
2003b). Thus it is not surprised to see the low mass galaxies 
have a different mass assembly mode. In fact, observations 
show that the evolution of low mass galaxies can be better in¬ 
terpreted in an "outside-in" framework. (Gallart et al. 2008; 
Zhang et al. 2012) 

It is now clear that the galaxy assembly mode has depen¬ 
dence on its already formed M*. Nevertheless, we know less 
about whether there is/or not a mass threshold which can 
roughly separate these two distinct modes. Answering this 
question requires the investigation of a sufficient large sam¬ 
ple, which expanding a wide mass range. Recently, using 
the lES data of 105 galaxies observed by the CALIEA project 
(Sanchez et al. 2012), Perez et al. (2013) found that galax¬ 
ies that are more massive than 5xlO° °M 0 grow inside-out, 
while lower mass galaxies grow outside-in. However, this 
sample is still relatively modest, and there are only about 20 
galaxies with M* < 10'°Mq galaxies in this sample (see their 
Figure. 1). To study a larger sample, in this Letter, we use the 
GALEX (Martin et al. 2005) photometry, combined with the 
SDSS (York et al. 2000) data, to investigate how galaxy as¬ 
sembly mode depends on stellar mass. Throughout this Letter, 
we assume a concordance ACDM cosmology with Um = 0.3, 
Ha = 0.7, Hq = 70 km s“' Mpc“', and a Kroupa (2001) IMF. 

2. METHOD AND DATA USED 

In this Letter, we will investigate the stellar age in the cen¬ 
tral region, as well that of the global, in individual galaxies. 
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For a galaxy, by comparing these two quantities, we can inves¬ 
tigate whether its central part is first assembled. Traditionally, 
stellar age is determined using broad-band optical photometry 
when spectral information is not available. However, the use 
of optical photometry has been plagued by the well-known 
age-metallicity degeneracy (AMD) whereby young, metal- 
rich stellar populations produce optical colors which are in¬ 
distinguishable from those produced by old, metal-poor pop¬ 
ulations (Worthey 1994). To avoid the AMD, We use the 
NUV-r color index as an stellar age indicator, since it is sen¬ 
sitive to low level recent star formation (RSF), but is unaf¬ 
fect by the AMD when stellar age is < 1 Gyr (Kaviraj et al. 
2007a,b). 

Following Wang et al. (2010)’s pipeline, we have con¬ 
structed a UV-optical matched photometric catalog. This 
catalog contains about 220,000 galaxies with uniform pho¬ 
tometric measurements on the resolution and PSF-matched 
GALEX+SDSS images, which are cross-matched by the 
SDSS DR8 (Aihara et al. 2011) spectroscopic galaxies and 
the GALEX GR6 database. For each galaxy, the fluxes were 
measured over 5 different apertures, with r =[1.5, 3.0, 6.0, 
9.0, 12.0] ", in the Ft/y,Af/y,M,.g,r,!,z bands. Along with 
the aperture photometry, we also measure the total magni¬ 
tudes of the galaxies with SEXTRACTOR (Bertin & Arnouts 
1996). All the magnitudes have been corrected for galactic 
extinction using the galactic dust map (Schlegel, Finkbeiner 
&L Davis 1998). The data reduction procedure is refereed to 
Wang et al. (2010) and Pan et al. (2014). 

To measure the NUV-r color index in the central region of a 
galaxy, one would expect the galaxy to be spatially resolved in 
both the GALEX and SDSS images. The GALEX NUV image 
has a pixel size of 1 pixel=l".5 and a point spread function 
(PSF) with full width at half-maximum (FWHM)=5".3. The 
pixel size and PSF of the SDSS image are Q"396 and 1".4, 
respectively. To ensure that the measured central NUV-r 
color index is not significantly affected by the relative poor 
resolution of the GALEX images, the chosen central aper¬ 
ture should not be smaller than the NUV PSF, e.g, rcentrai 
should not be smaller than 5"3/2 « 2".7. In the 5 apertures, 
the r = 3".0 aperture is thus most suitable to be chosen as 
the "central aperture". Then, we select the galaxies passing 
the following as our parent sample: 1) SDSS minor-major 
axis ration b/a > 0.5 ;2) SDSS Rgo > 5.0", where Rgo is 
the radii enclosed 90% of the SDSS r-band petrosian flux; 

3) z = [0.005,0.05], where z is SDSS spectroscopic redshift; 

4) mNuv < where niiimit = 23.0 mag is the limited mag¬ 

nitude in the NUV band. The limit on b/a is to minimize 
the dust reddening effect on our results. Finally, we limit 
these galaxies to have stellar mass M* > 10®'° Mq, to form 
a volume-completed sample at z < 0.05 (Schawinski et al. 
2010). The final sample size is Vgai=ll, 294. The stellar 
mass M*, and the spectral index D„4000 used in this Letter 
was drawn from the JHU/MPA database 

In the upper pannel of Figure. 1, we show the central flux 
fraction (/central) distribution for our sample, where /central was 
defined as /central = Flux^^^n FluxiotuX- R can be seen that in 
the SDSS r-band, the median value of /central is ^ 0.25. One 
will also find that the /central distribution is much broader in 
the NUV band. The Kolmogorov-Smirnoff (KS) test on these 
two distributions gives a small probability (<0.01) that they 
are drawn from the same distribution. This may be due to 

* http://www.mpa-garching.mpg.de/SDSS/DR7 
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Fig . 1.— Upper panel: The /central distribution for our sample, where /central 
was defined as /central = Flux^_y/ ^/FIuxiova- The GALEX NUV and SDSS 
r band is shown in blue and red histograms, respectively. Lower panel: the 
cumulative curve of /central- 

the highly different RSF in the galaxy central regions. Fig¬ 
ure. 1 illustrates that, even for the Rgo ~ 5.0" galaxies, in the 
SDSS r-band, their /central not higher than 0.4. Thus for 
this galaxy sample, our chosen central aperture should be ap¬ 
propriate. 

3. RESULTS 

3.1. The color—magnitude diagram 

The color-magnitude diagram, or the color-mass diagram 
(CMD), is widespread used in the literature to diagnose the 
galaxy evolutional stage (Faber et al. 2007; Kaviraj et al. 
2007b; Schawinski et al. 2014). Based on their positions 
on the CMD, galaxies can be roughly categorized into two 
main populations, e.g., the red sequence galaxies and blue 
cloud galaxies. Galaxies that lie between the "blue cloud" 
and "red sequence", are called the "green valley " (GV) galax¬ 
ies (Strateva et al. 2001; Baldry et al. 2004). GV galaxies 
were traditionally thought to be the transition populations be¬ 
tween the star-forming and the quenched galaxies (Pan et al. 
2013; Mok et al. 2013; Schawinski et al. 2014; Pan et al. 
2014; Vulcani et al. 2015). To show the CMD clearly, we 
have divided the sample into 3 subsamples, according to their 
A NUV-r, where ANUV-r = (NUV-r^j,, Q)-(NUV-rtotai)- 
NUV-r^_^// and NUV-rtotai are the color indices measured 
in the central r = 3".0 aperture and that of the whole galaxy, 
respectively. The typical error of A NUV-r, ctq, is derived as 
Jq" N(a)/Ngsxdcr=0.S0, where N(a) is the number of galax¬ 
ies with a A NUV-r error of a. For this sample, cto is 0.15 
mag. In what follows, we called galaxies with A NUV-r>2(To 
as "red-cored" galaxies. The "blue-cored" galaxies are those 
with A NUV-r<-2tTo- They are denoted in red and blue sym¬ 
bols in panel a) of Figure. 2, respectively. Those galaxies 
with flat color gradients, e.g., galaxies withIA NUV-rl<2tTo, 
are shown in green symbols. All AGNs in our sample are 
narrow-line (obscured) AGNs, so there should be no contri¬ 
bution from AGN continuum to the NUV-r color index. 

In panel b), panel c) and panel d) of Figure. 2, we show 
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Fig. 2.— a):The relation between NUV—r _j// ^ and NUV—notal, where NUV-r^_j// ^ and NUV-rtotal are the color indices measured in the central r = 3”.0 
aperture and that of the whole galaxy, respectively. ANUV—r is defined as ANUV—r = (NUV-r^ j// q)—(NUV— ftotal)- b): The M* .vs NUV-r (total) diagram for 
galaxies with ANUV-r < —0.3. c): CMD for galaxies with -0.3 < A NUV-r < 0.3; d): CMD for galaxies with A NUV-r > 0.3. The two pink lines denote 
our "green valley" definition. 



Fig. 3.— The fractions of blue-cored (blue symbols), flat-gradient (green symbols) and red-cored (red symbols) galaxies as a function of M*. From left to 
right, we show the results for the blue cloud, green valley and red-sequence. Errors are derived from 100 bootstrap resamplings. 


the M*- (NUV-r) diagram for the blue-cored, flat-gradient 
and red-cored galaxies, respectively. The gray scale repre¬ 
sents the galaxy number density of the whole sample. On 
the CMD, one can see the evident "color bimodality". The 
GV region is arbitrarily defined by the two pink lines, with 
0.6 X M* -2.4<(NUV-r)Gv< 0.6 x M* - 1.4, where M* is 
logarithm stellar mass. As shown in Kaviraj et al. (2007a), 
galaxies with an NUV-r color less than 5.0 are very likely 
to have experienced recent star formation. In the blue cloud 


and GV region, the "blue-cored" phenomenon is thus corre¬ 
sponding to centrally concentrated star formation. In panel c), 
we show that flat-gradient galaxies distribute similarly with 
the whole sample on the CMD. This is expected since they 
are most numerous. Strikingly, panel d) shows that there 
are nearly no "red-cored" galaxies in the GV region below 
M* - 1O*°°M0. 

In Figure.3, we show the fractions of the 3 kind galaxies, 
scaled with M*, in the blue cloud, green valley and red se- 
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Fig. 4. — panel a):The central Dn4000 .vs NUV-r (total) relation for the whole sample, panel b), panel c) and panel d) show this relation for the log 
M* /Mq<10.0, log M* /Mq=[10.0,10.5] and log M*/Mq>10.5 subsamples, respectively. Color symbols are the same in Figure. 2. 


quence. It is clear from Figure. 3 that at > 
the blue cloud and GV region is mostly populated by red- 
cored galaxies. In the GV region, we see a significantly 
increase of blue-cored M* < 10 '*’ galaxies when com¬ 
pared to the blue cloud, at the expense of a decrease of red- 
cored and flat-gradient galaxies. This can be interpreted as 
the outside-in quenching processes (regardless of its working 
mechanisms) transform a significant fraction of the red-cored 
and flat-gradient galaxies into blue-cored ones, and drive them 
evolve form blue cloud to GV. In the red sequence, galaxies 
are mostly with flat color gradients. This is expected since 
when NUV-r color index is no longer sensitive to stellar age 
older than 1 Gyr (Kaviraj et al. 2007b). To summary. Figure. 
2 and Figure. 3 confirms the findings of Perez et al. (2013) 
that < 10*° galaxies mainly grow "outside-in". The 
importance of the "inside-out" mode steadily increases with 
M*, and becomes dominant at around M* > lO'^'^M©. 

3.2. The D„4000 versus NUV—r diagram 

In this part, we will use a spectroscopic index, the D„4000 
index measured on the SDSS r = 1".5 aperture spectra as an 
central stellar age indicator. The Z)„4000 is widely used as 
a mean stellar age indicator in the literature, since it is sen¬ 
sitive to the ratio of present- to past- average star formation 
rate (Kauffmann et al. 2003a; Brinchmann 2004), but is less 
affected by dust and metallicity. 

In panel a) of Figure.4, we show the D„4000 versus NUV-r 
(total) relation for our sample. The grayscale represents the 
number density of galaxies. We can see that there is a positive 
relation between D„4000 and NUV-r (total), as also found in 
Kauffmann et al. (2007). The scatter of this relation reaches 


its maximum at NUV-r « [3.0,5.0], .e.g. the GV region of 
the CMD. In panel b), panel c) and panel d), we show the 
the D„4000 versus NUV-r (total) relation for three subsam¬ 
ples, which are binned according to their M*. The color coded 
symbols are the same as those in Figure. 2. In Figure. 4, one 
can find that the red and the blue symbols are still separable. 
This confirms that the NUV-r _j// is equivalent to the SDSS 
D„4000 index. 

In panel b) of Figure. 4, we find that, nearly all of the 
M, < 1 O*°°M 0 galaxies have moved to the red sequence 
when their central D„4000>1.6. In contrast. In panel d), we 
find that a large portion of the M* > 1 O*°'^M 0 galaxies still 
lie on the blue cloud or the GV region when their central 
D„4000>1.6. Panel b) and panel d) together show that, there 
seems to be two separable evolution paths through which the 
blue cloud galaxies evolve into the red sequence. In this work, 
it is naturally to interpret these two paths as the "outside-in" 
and the "inside-out" path, respectively. 

In panel b), we show that the flat gradient galaxies are on 
the same path as those of "blue-cored" galaxies. Thanks to 
the the accurate measurements of D„4000 and NUV-r (total), 
we can now conclude that at M* < 1 O*° °M 0 , the flat gradi¬ 
ent galaxies follow a similar evolutionary path of the "blue- 
cored" when they move to the red sequence, e.g., they evolve 
"outside-in". 

A potential issue which may affect the result of Figure. 4 
is that we haven’t applied dust reddening correction to the 
NUV-r (total) color index. Since the color index of the 
inclined galaxies are more easily reddened by the galactic 
dust lanes, to assess the role of dust on our result, we show 
the D„4000 versus NUV-r (total) relation for the fi/a < 0.5 
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Fig. 5.— Same to Figure. 4, but shown for the inclined galaxies. 

galaxies in Figure. 5. by comparing Figure. 4 and Figure. 5, 
we find that the result is not changed for the M* < 10 “’°Mq 
galaxies. At D„4000 > 1.6, many M* > 10*°'^ Mq galaxies 
are reddened by dust, but their distributions can still be sepa¬ 
rated from those of the M* < 10*° on this diagram. This 
test suggests that the two separable paths revealed in Figure. 
4 is intrinsic. 

4. SUMMARY AND DISCUSSION 

In this Letter, we investigate how galaxy mass assembly 
mode depends on stellar mass M*, using a large sample of 
'^lO, 000 low redshift galaxies. The use of GALEX photom¬ 
etry and SDSS data enables us to study a large sample with¬ 
out the IFS data. We compare the stellar age indicator mea¬ 
sured in the galactic centers to that of the whole galaxy, to 
determine whether the galaxy grows "outside-in" or "inside- 
out". We find that: in the M* -( NUV-r) green valley, the 
M* < 10'° Mq galaxies are mostly "blue-cored" or have flat 
color gradients, whereas the M* > 10'°'^ Mq galaxies are 
"red-cored". When their central D„4000 index values exceed 
1.6, the M* < 10'° ° Mq galaxies have lied on the UV red se¬ 
quence, whereas a large fraction of the M* > 10*°'^ Mq galax¬ 
ies still lie on the UV blue cloud or the green valley region. 
These findings suggest that the main galaxy assembly mode 
is transiting from "the outside-in" mode to "the inside-out" 
mode at M* < 10'° Mq and M* > 10'°-^ Mq. 

We interpret these results as a compromise between the in¬ 
ternal and external process that driving star formation quench¬ 
ing in galaxies. The internal process, like the AGN feedback 
(e.g., Croton et al. 2006) and the supernova feedback induced 
from the strong starburst activity (Geach et al. 2014), can blow 
the gas out of the galactic center, causing the "inside-out" evo- 


D4000 (r<1.5") 

lution. The external process, like the environment effects, 
will strip the gas content of galaxies firstly from their out¬ 
skirts, causing the "outside-in" evolution. The capability that 
a galaxy can retrieve the lost gas, or accreting new gas that 
rebuilds its outskirts, is its gravity. It is thus not surprising 
to see the galaxy assembly modes have strong dependence on 
M*. 

Panel b) of Figure. 2 shows that a significant fraction of the 
M* < 10*°^ Mq galaxies are "blue-cored". Previous works 
found that galaxies of this type are usually bulge-dominated 
systems, like the Eh-A galaxies (Yang et al. 2008) or the early 
type galaxies (Lisker et al. 2007; Suh et al. 2010; Pan et al. 
2014). This phenomenon can be either explained as the "sur¬ 
vived" star formation of an late type galaxy which is suffering 
gas stripping, or explained as the "excess" central star forma¬ 
tion, probably triggered in an merge event. Since the merge 
rate is rare in the local Universe, we suggest that the "survived 
star formation" explanation is the case for the majority of the 
"blue-cored" galaxies, especially for the low mass ones. 

M* < 10'° Mq galaxies are mostly satellites of massive ha¬ 
los. In a group/cluster environment, ram pressure stripping 
(Gunn & Gott 1972) and tidal interactions, can remove the 
gas from galaxies, especially for those with low stellar mass 
densities (Zhang et al. 2013). Other gentler processes such as 
"Starvation" and "Harassment", are also at work (Barazza et 
al. 2002; Weinmann et al. 2009). The important role of dense 
environment on the formation of red dwarf galaxies is con¬ 
firmed by the works focused on the Virgo and Coma clusters 
(see Boselli & Gavazzi (2014) for a review). In lower halo 
mass, Geha et al. (2012) found that quenched M* < 10° Mq 
galaxies do not exist in the field, suggesting that the internal 
process alone is not sufficient in quenching star formation of 
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low mass galaxies. At higher redshifts, using the COSMOS 
data. Pan et al. (2013) found that the important role of en¬ 
vironment on the evolution of M* < 10'° Mq galaxies have 
been in place at z < 0.7. 

The hierarchical paradigm for galaxy formation implies that 
galaxies experience a variety of environments during their 
evolution, and may have complex SFH at their early assem¬ 
bly epoch. In Figure. 2 and Figure.4, we can find that there 
are many "red-cored" dwarf galaxies in the blue cloud. How¬ 
ever, these galaxies are not quenched, even their cores are 
older than the global. It is worth to mention why we focus 
on the GV galaxies at this point. GV galaxies are more close 
to the "fully assembled" stage than the blue cloud galaxies 
(here we do not include those GV galaxies due to "rejuve¬ 
nation"), hence holding important clues to the galaxy final 
assembly process. As can be seen in panel a) of Figure. 4, 
when galaxies evolved into the GV region, there seems to be 
two separable paths connecting the blue cloud and the red se¬ 
quence. This is expected if there is a dichotomy in the galaxy 
assembly modes. 

We find that the log M*/M0=[1O.O,1O.5] regime is a tran¬ 
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